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Abstract 

A general method to construct basis functions for fcrmionic systems which ac- 
count for the SU(2) symmetry and for the translational invariance of the Hamil- 
tonian is presented. The method does not depend on the dimensionality of the 
system and it appears as a natural generalization of the Bethe Ansatz to the case 
of non integrable systems. As an example we present the block diagonalization 
of the Hubbard hamiltonian for finite number of sites in one and two dimensions. 
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During the past year a great deal of interest has been devoted to strongly correlated 
Fermi systems in two spatial dimensions because of their possible role as models for high T c 
superconductivity pj. Among these systems, the Hubbard model is certainly the simplest 
non trivial model for interacting electrons in a solid. Its Hamiltonian consists of a kinetic 
part, representing the electronic hopping between next neighbor sites and a a part represent- 
ing the repulsive Coulumbian interaction of electrons with opposite spins on the same site. 
In spite of this apparent simplicity, the mathematical and the physical properties of this 
model in dimensions higher than one (magnetic properties of the ground state, existence of 
possible Mott transitions between insulating and conducting states, off-diagonal long range 
order, etc.), are still poorly understood. On the other hand, the large dimension of the 
Hilbert space of the states makes any numerical attempt to compute the spectrum of these 
systems possible only for clusters of very small size. A cluster of 6 x 6 lattice sites for the 
Hubbard model is already behind the computational power (using algorithms which account 
for the translational symmetry and for the conservation of the z component of the total 
spin S of the system), of any supercomputer nowadays available. The main difficulty is the 
lacking of general methods which allow to account for all the symmetries, both continuous 
and discrete, of these systems. For bosonic systems, basis functions which account for the 
conservation of the number operator and for the translational invariance on a periodic lattice 
can be constructed by the so called number state method 0] . In the case of the 2D Hubbard 
model one would like to use as basis functions the simultaneous eigenfunctions of S" 2 , of the 
total number electrons N and of the translational operator on the lattice. With respect to 
these functions the Hamiltonian would acquire a block diagonal form with blocks of minimal 
size (if all the symmetry is included), reducing in a significant manner the memory needed 
for the diagonalization procedure. 

The aim of this letter is to present a general method to construct basis functions for 
strongly correlated Fermi systems which span the irreducible representations (irreps) of 
SU(2), U(l), and of the translational group on a periodic lattice. The method does not de- 
pend on the dimensionality of the system and it appears to be a natural generalization of the 
Bethe Ansatz to the case of non integrable systems (in the case of the one-dimensional an- 
tiferromagnctic Hcisenberg chain it is completely equivalent to Bethe Ansatz). The method 
is based on the following points: i) The invariance of the total spin of the system under the 
permutation group Sf allows to construct basis functions which span the irreps of SU{2). 
This is done by using "filled" Young tableaux for fermionic systems recently introduced in 
rcf. H. ii) Any discrete group is a subgroup of the permutation group, thus one can project 
the above eigenfunctions on the Sf subgroup of interest (in the present case on the sub- 
group corresponding to the translations on the lattice). The commutation of the translation 
operator with the total spin S, assures that the projected functions will be simultaneous 
eigenfunctions of both operators. We illustrate the method by taking as working example 
the Hubbard model both in one and in two dimensions. 

The Hubbard Hamiltonian is written as 



As usual, U represents the onsite Coulomb repulsion energy between electrons, and n il 
are respectively, spin up and spin down occupation number operators at site i, t denotes the 
hopping energy between sites, a denotes the possible spin states of the electrons (i.e. a =f 




(1) 
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or J,), and cj , Cj a are fermionic creation and annihilation operators 

{ci^Cjv} = {cti„> c L |=0, |c v , c] ff , | = SijS^. (2) 

As well known the Hamiltonian ([]]) posses, besides the translational symmetry and the SU(2) 

invariance under rotations of the total spin S of the system, also the U(l) conservation of 
the total number of electrons 

^EK+^i)' ( 3 ) 

3 

and the so called particle-hole symmetry. The conservation of N and S z obviously implies 
the separate conservation of the number of electrons with spin up = ■ n JT and of 
electrons with spin down N± = J2j n ji ■ To define the Hilbert space of the states we denote 
with |3 >, |2 >, |1 >, |0 > the possible states on a generic site i.e., respectively, the doubly 
occupied state, the single occupied state with spin up, the single occupied state with spin 
down and the vacuum (the numbers 3,2,1,0 are here used as quantum numbers to characterize 
them). Since there are four possible states at each site, the dimension of the Hilbert space 
for a lattice with / sites is just 4/. On the other hand the conservation of the number 
operator N allows to decompose the Hilbert space K into a direct sum of eigenspaces Kn 
corresponding to a fixed value of N. The dimension of these spaces is just the number 
of ways N electrons can be placed in 2/ boxes i.e. djy = pf'-Ny.N'. • ^ ms ^ e& ds to the 
diagonalization (for each N) of a finite d,N X c/jv matrix (we can restrict N to N < f since 
/ < N < 2f follows from the particle-hole symmetry). We then use the invariance of S 2 
under Sf to construct S 2 ,N eigenfunctions which allow to further block diagonalize each 
djv x d-N matrix with respect to the irreps of SU (2). To this end we recall that the irreducible 
representations of Sf can be labelled by all possible partitions [fi-.-fk] of / into k parts, 
with fi integers obeying /i > / 2 > ... > fk and /i + /2 + — + /fe = /■ To each partition it is 
associated a Young tableau of type {/i, f%, ...} i.e. having f\ boxes in the first row, ji boxes 
in the second row, etc., each tableau corresponding to different symmetry class operations 
oiSfg. 

In order to obtain the highest weight vectors of SU(2) (i.e. eigenvectors of S 2 , S z 
belonging to S = S z ) with a given Sf symmetry we use the idea of Young tableaux "filled" 
with quanta, introduced for bosonic systems in ref. j^] and extended to the fermionic case in 
ref . H . To this end we observe that for a fixed TV = N-\ + N± , the possible values of iVf , N± 
compatible with it are all the partitions of N into two parts, each partition being associated 
with a well defined value of S z . Let us introduce the quantum number 

M = 3N l + 2(N 1 ~N l ) (4) 

and consider all the partitions (mi,m2,...,m/) of M into / parts with m, = 0,1,2,3 (in 
eq. (^) and in the following, we restrict to iVf > iVj, since the cases N-\ < N± follows from 
these by interchanging spin up with spin down electrons). Since the order of the quantum 
numbers m; in the partition is unimportant we fix it to be mi > m,2 > ... > mf. We remark 
that each partition (mi,Tf>2, ...,m/) of M is associated with an eigenstate |mi,rri2, > 
of S z . This leads to a family of states organized into levels. From each M level we construct 
eigenmanifolds of S 2 with a given Sf symmetry by filling the quanta n%i, characterizing that 
level in the boxes of a Young tableau according to the following rules: i) The quanta must 
be not increasing when moving from left to right in each row or when moving down each 
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column of a given tableau, ii) The quanta referring to spin up and spin down states (i.e. 
mi = 1,2) cannot appear more than once in a row. iii) The quanta referring to doubly 
occupied states or to empty states (i.e. m, = 3, ) cannot appear more than once in a 
column. These rules directly follow from the permutational properties of the four states 
1 3 >, 1 2 >, |1 >,|0 > and from the symmetrization and antisymmetrization property of, 
respectively, rows and columns of a given Young tableau. By using these rules we construct, 
for each M level, a family of filled Young tableaux for each Sf symmetry. To pass from filled 
tableaux to states we apply Young symmetrizer and antisymmetrizer operators which take 
into account the Pauli exclusion principle, i.e. every time the symmetry operator involves 
the permutation of n a spin up electrons and spin down electrons an extra (— l) n »+™*> 
factor is included (note that due to the permutational properties of the \2 > and |1 > states, 
one must count only the interchanges of spin electrons and spin down electrons separately 
involved in each permutation). These states, by construction, are eigenstates of S z with a 
given Sf symmetry but, in general, they are not eigenstates of S 2 . In order to identify the M 
levels with the eigenmanifolds of S 2 we must characterize the filled tableaux corresponding 
to highest vectors of SU(2) (highest weight filled tableaux). This can be done by noting 
that a change of a 1 into a 2 (1-2 flip) in a filled tableau corresponds to increases M by 1 
i.e. to pass to a filled tableau of the M + 1 level. We have therefore, that the filled tableaux 
which survive 1-2 flips (i.e. the one that satisfy the filling rules also after a 1-2 flip) are 
the ones for which S > S z . In this way one "extract" from each M level the highest weight 
vectors of SU(2). We remark, however, that by a 1 — 2 flip two different M tableaux may 
be associated with the same M + 1 tableau. In this case one can prove that the linear 
combinations tp± — tpi ± ip2 of the states corresponding to the M filled tableaux, produce 
one S = S z state ((f) +) and one S = S z + 1 state ((f)-)- Note that these states are not 
necessarily orthogonal so that, in general, a final Gram-Schmidt orthonormalizcd procedure 
must be applied. We remark that the functions so constructed are good basis functions to 
solve fermionic systems with infinite-range interactions (in the thermodynamic limit this 
should correspond to exact mean- field calculations). The application of these functions to 
the Hubbard model with unconstrained hopping is discussed in rcf. (3). By using these 
functions one can easily get the following characterization of the ground state for the Sf- 
invariant Hubbard system for t > 0: the N=l,2 ground state is always associated with a 
tableau of type {/} while for / > 3 it is associated with tableau of type {/— (N— 2), 2, ...,2} 
for N even or of type {/ — (N — 2), 2, 2, 1} for N odd. Furthermore, in the ground state 
S has always its minimal value i.e. S = for N even or S = i for N odd. 

We now come to the problem of projecting the above functions on the subgroup, say G, 
of Sf of physical interest (for translations on a ID periodic lattices G is just the abelian group 
Cf corresponding to the cyclic permutations). This problem is similar to the one encountered 
in perturbation theory when a perturbation reducing the symmetry induces a splitting in 
the energy levels. Let us denote by D(R),R € Sf the irreps of Sf. A representation of G 
is readily obtained by selecting among the matrices D(R) those corresponding to elements 
of G. These representations however arc in general reducible i.e. they can be expressed in 
terms of irreps D'"' of G as D(R) = ^ v c v D^ v \R) with c v non negative integers counting 
the number of times D^ v ' appears in D. By denoting with gi the number of elements in the 
class Ki of G and with g the order of this group, one easily express the integers c u in terms 
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of the characters Xi X °f respectively Sf and G as 

<V = ~^T f 9iXi' ) *Xi- (5) 

" i 

This gives the splitting of the irreps of Sf (i.e. of the Young tableaux of a given type) into 
the irreps of G (star in Eq.(|J) denotes complex conjugation). The eigenfunctions ip of N, S 2 , 
corresponding to the above highest weight filled Young tableaux, are projected on the v— th 
irrep of G by using the projection operator defined by 

ipM = P (v H = — x (v) *{R) U R ■ ip (6) 
g R 

where the sum is over all the elements R of the 5/ subgroup, n v is the dimension of the i^-th 
irrep of G, (R) the corresponding characters and Ur the operator associated to the group 
element R Q . By taking G to be the subgroup corresponding to the lattice translations T n 
and by projecting all the functions corresponding to highest weight filled Young tableaux 
with a given value of S and N, we get the simultaneous eigenfunctions of N,S 2 ,T n with 
respect to which the Hamiltonian aquires block diagonal form. 

Let us illustrate the method with an explicit calculation on the Hubbard model with / = 4. 
We consider the four sites disposed in two different configurations: the first corresponding 
to a 1-D periodic chain, the second to a 2-D periodic square lattice. In the first case the Sf 
subgroup of interest is the cyclic group C4, while in the second case is the group C2/1. They 
are both abelian groups with one dimensional irreps. Let us denote with A, B, E\, E2, the 
irrep of C4 and with A u , A g , B u , B g the irrep of C%h (we refer to the standard notation of 
point-symmetry groups) . The two irrep E\ , E% of the group C4 are one the complex conjugate 
of the other thus they physically correspond to a double degenerate level (this is true also 
for other C / groups) . This accidental degeneracy is connected to time reversal invariance of 
the Hubbard Hamiltonian i.e. the complex conjugate of an eigenfunction is automatically 
an eigenfunction with the same energy. We take advantage of this fact by considering E\ , Ei 
equivalent to a single irrep E of dimension two. For brevity we concentrate only on the case 
S = at half filling (N = 4) (a detailed analysis will be published elsewhere Q) . In Table 1 
we report all the S = highest weight filled tableaux together with their splittings in terms 
of the irreps of the groups C4, and C^h- From this table we see that for the 1-D chain one 
gets two blocks of dimension 6x6 associated with the irrep A and B (this giving 12 S = 
nondegenerate eigenvalues), and one 4x4 block associated to the irrep E (giving four doubly 
degenerate eigenvalues). In the case of the 2D chain we see that the accidental degeneracy 
in the E representation is removed, and we have three 4x4 blocks (respectively associated 
to the irreps A u , B u , B g ) and one 8x8 block associated to the A g representation. Let us 
concentrate here only on the ground state. To this end we remark that for the S'j-invariant 
Hubbard system the ground state is characterized by a S = Young Tableau of type {2, 2} 
(see above discussion and ref. || for details). We conjecture that the projection on the 
translational subgroup will not alter this situation, i.e. that the ground state belongs to 
one of the irrep in which the Sf ground state splits (we think this conjecture holds true for 
general cases if / is even and t > 0) . This is indeed what happens for the present cases. We 
find that the ground state is of type B for the 1-D chain and of type A g for the 2-D chain. 
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The block associated to the B representation of the ID chain is given by 
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while the block associated to the A g representation of the 2D chain is: 
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To check these results we have numerically diagonalized H in both cases in the N = 4 
eigenspace of dimension — 70. From these calculations it follows that, for t = 1, U = 2 and 
for f = N = 4, the ground state of the ID chain is non degenerate with energy E — —2.82843 
while for the 2D chain the ground state is non degenerate with energy E = —6.681695. One 
easily verify that these values coincide with those obtained by diagonalizing the blocks 
respectively in Eq. (|7j) and Eq.(||). 

In closing this letter we remark that the above method of block diagonalizing H is quite 
general and can be applied to more complicate fermionic systems such as the Anderson model 
as well as to other subgroups (besides translations) of physical interest such as the invariance 
group of fullerene molecules. Furthermore we note that the study of the Heisenberg model 
and of the t-J model in one and two dimensions directly follow from the above analysis 
by restricting the space of the single site states respectively to |0 >, |1 > (Heisenberg) or 
|0 >, |1 >, 1 2 > (t-J). We also remark that the method is completely algebraic and can be 
easily implemented on a computer (by using Mathematica we have set up packages which 
performs all the operations to construct the above basis functions). We hope our method 
will facilitate future numerical studies of strongly correlated Fermi systems contributing to 
the understanding of the physical properties of these systems in two and three dimensions. 
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Table Captions 



Table 1 

Decomposition of the filled Young tableaux corresponding to S = highest weight vectors 
of SU(2) for / = N = 4, in terms of the irrep of the groups C4, and C^h- The sum of 
tableaux denotes the (plus) linear combination of the corresponding states. 
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